Roughness-induced critical phenomena in a turbulent flow 
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I present empirical evidence that turbulent flows are closely analogous to critical phenomena, from 
a reanalysis of friction factor measurements in rough pipes. The data collapse found here corresponds 
to Widom scaling near critical points, and implies that a full understanding of turbulence requires 
explicit accounting for boundary roughness. 
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Turbulent phenomena in fluids |lj are characterized by 
strong fluctuations and power-law spectraQ that are 
suggestive |3J of the power-law correlations observed in 
critical phenomena at continuous phase transitions Q. 
Attempts to link these two sets of phenomena have pri- 
marily, but not exclusively^ |fj, focused on the calcu- 
lation of scaling exponents |jj in ideal systems which are 
isotropic and homogeneous, neglecting boundaries. How- 
ever, no analogue has been found for the wide variety of 
thermodynamic scaling phenomena 8] that are an equally 
fundamental aspect of continuous phase transitions, and 
whose elucidation |9j led to a complete understanding of 
critical phenomenaflCj] . Here we consider the important 
role of boundary roughness on fluid flow, by reanalyzing 
Nikuradse's experimental measurements |f 1 \ of the fric- 
tion exerted on a turbulent fluid by the walls of a rough 
pipe. We show that the data as a function of Reynolds 
number and relative roughness collapse onto one univer- 
sal curve, when appropriately scaled. This analogue of 
Widom scaling 8] implies that boundary roughness must 
be included in a complete description of turbulence, and 
establishes the long sought-after precise connection to 
critical phenomena. 

Turbulent flows are characterised by their Reynolds 
number, defined as Re = UL/v, where U is a typical 
velocity at the length scale L, and u, the kinematic vis- 
cosity, is the viscosity of the fluid divided by its den- 
sity. In 1941, Kolmo gorov [l2] and Obukhovp^. recog- 
nised that at large enough Reynolds numbers, fluid mo- 
tion is, over a wide range of length scales, a dynamical, 
energy-conserving but irregular swirling motionQ] gov- 
erned by inertia, rather than a dissipative phenomenon. 
Thus, they pointed out that for this so-called inertial 
range of scales, observables should be independent of 
v. In particular, for the inertial range, the turbulent 
energy spectrum E of longitudinal velocity fluctuations 
5vk in wavenumber space, k, can only depend upon the 
mean energy transfer rate e and k itself in a manner dic- 
tated by dimensional analysis: E(k) = (\6vk\ 2 ) = efc~ 5 / 3 . 
This experimentally-verified^^ power-law scaling (often 
referred to as K41) applies on small scales in a turbu- 
lent flow, but not so small that molecular viscosity be- 
comes important. The existence of a wide range of length 



scales, over which power-law (and thus scale-invariant) 
correlated fluctuations are found, is reminiscent [2} of the 
power-law fluctuations on many length scales that ac- 
company critical phenomena 4]: for example, in a fer- 
romagnet near its critical point, the Fourier compo- 
nent of the magnetisation M at wavenumber k satisfies 
G(k) = (|Affc| 2 ) cx k~ 2+v , where n is an anomalous scal- 
ing exponent that describes departures from mean field 
theory. This power-law scaling applies at large scales, 
and is independent of the small scale details of the sys- 
tem, such as the nature of the crystal lattice. 

Power-law scaling of correlation functions is, however, 
only one of two key aspects of critical phcnomena[4j . The 
other, equally important aspect is the phenomenon of 
data collapse, or Widom scaling 8J: for example, in a 
ferromagnet, the equation of state, nominally a function 
of two variables, is expressible in terms of a single re- 
duced variable that depends on a combination of external 
field and temperature. What is the turbulent analogue 
to Widom scaling in thermodynamics? To address this 
question, it is necessary to examine data on the large- 
scale properties of turbulence, for example the friction 
factor in pipe flow. 

In 1932 and 1933, Nikuradse undertook a seminal se- 
ries of measurements of flow in nominally smooth [l^ and 
rough pipes [lj, measuring inter alia the friction factor 
/, related to the pressure drop across the piped 113. 
These measurements have remained the benchmark in 
the field, being the only systematic measurements of a 
single flow geometry over such a wide range of Reynolds 
numbers. Nikuradse's experiments on sand-roughened 
pipes used sand grains of a well-defined size r, repeated 
over a wide range of values of r and pipes of different 
radii D. Nikuradse was able to verify the expectation of 
hydrodynamic similarity: the flow properties depend on 
the roughness only through the combination r/D. Niku- 
radse presented results for the shear force per unit area 
r exerted by the flow on the walls of the pipe in the form 
/ = r/pU 2 , and these data are plotted in figure 0). 

There are several broad features to note in these data. 
First, in the smoothest pipes studied, for Re < 100,000, 
in the turbulent regime, the friction factor is a decreas- 
ing function of Re, varying to a good approximation 
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FIG. 1: (Color online) Friction factor for turbulent flow in 
a rough pipe, as reported by Nikuradse 11]. The data were 
measured at different values of Re and r/D, and extracted 
from Nikuradse's tabular and graphical presentation. 

in a manner usually attributed to Blasius|0, 0] as 
/ - Re- 1 / 4 . As the roughness scale is increased, and 
thus for smaller values of D/r, the extent over which the 
Blasius scaling extends becomes smaller. Thus, we can 
represent this feature by the statement that / ~ Re^ 1 / 4 
asymptotically as r/D — > 0. Asymptotically for rough 
pipes at large Re, the friction becomes independent of 
Re, and depends only on the roughness r/D, varying to 
a good approximation, according to Strickler's law|20j. as 
/ ~ (r/D) 1 / 3 . These broad characteristics are also visi- 
ble in aggregate in other pipe flow datalHEllHElE^. 
but no other single data set captures the full range shown 
here, with as little scatter evident in the data. 

These features place strong constraints on the func- 
tional form of the friction factor f (Re, r/D). In fact, 
these constraints precisely parallel those on thermody- 
namic properties of ferromagnets near the critical tem- 
perature T c , as a function of reduced temperature t = 
\T — T c \/T c and external magnetic field H. For exam- 
ple, at the critical temperature, the magnetic equation 
of state has the form M ~ H 1 / 6 for t = 0, where 5 is a 
critical exponent whose value can be computed by renor- 
malisation group (RG) theory 0; and for zero field, the 
magnetisation continuously approaches zero as T — > T c ~ 
with a power law variation M ~ ir for H = 0, where ft is 
another critical exponent whose value can be computed 
by RG. Widomp discovered that these properties all fol- 
lowed if the thermodynamic free energy obeyed certain 
scaling properties, later shown to follow from renormal- 
ization group considerations |il llfJ|. 

We connect scaling in turbulence with that in critical 
phenomena by observing that the limit Re — > oo is anal- 
ogous to the limit t^ 1 — > oo, whereas the limit H — > 
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FIG. 2: (Color online) Friction factor for turbulent flow in a 
rough pipe, as reported by Nikuradse|llf . scaled according to 
the text. Inset: the same plot on a logarithmic scale. 

is analogous to the limit r/D — > 0. Thus, the analytic 
properties of the friction factor can be derived if we follow 
Widom's scaling argument Q, and propose a scaling form 
for the friction factor: f (Re, r/D) = Re- 1 l 4 g(Re a r/D), 
where g(z) is an unknown scaling function of a single 
variable z, which tends to a constant for small values of 
z, and a is an exponent that we can determine by the 
requirement that the Re dependence should cancel out 
of the formula for / at large Re, leaving the Strickler law 
/ ~ (r/D) 1 / 3 . This requires that g(z) ~ z 1 / 3 as z — > oo, 
and therefore a/3 = 1/4. Thus, we conclude that 

/ = Re- 1 l 4 g(Re 3 ' 4 (r/D)). (1) 

The scaling form of Eq. (JTJ predicts that the turbulent 
friction factor data, measured as a function of both r/D 
and Re, and thus in principle occupying a two dimen- 
sional space, will actually collapse onto a one-dimensional 
curve, when plotted as fRe ' versus (r/D)Re 3 / 4 . The 
test of this prediction is shown in figure , where Niku- 
radse's data occupying the plane of figure 0) collapses 
onto a single curve when plotted in the reduced variables 
of Eq. (JJJ . Note that the data collapse occurs for those 
data that lie between the Blasius and Strickler regimes 
only. Small deviations from the data collapse are visible, 
but it is not clear to what extent these reflect uncertain- 
ties in the data[2^| or something more fundamental. 

The scaling function that we have extracted from the 
data is unlikely to be universal in the sense of being in- 
dependent of the nature of the roughness of the pipe. 
By analogy with the effects of long-range interactions in 
critical spin systems 0, it seems probable that self-affine 
roughness [2^ will have a different effect on the flow than 
periodic single-scale roughness, and this can be reflected 
in the scaling function, the scaling exponents, or both. 
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Inspection of the data from different pipe flow experi- 
ments suggests that the first possibility is the most likely, 
but this remains to be checked in detail. 

The interpretation of the analogy discovered here is 
very natural. In magnetic systems, for example, the 
power-law scaling of G{k) is now understood to be as- 
sociated with the fact that magnetisation can be induced 
in a magnet by the application of an external magnetic 
field H . The sensitivity of the order in a magnet to a per- 
turbation by H becomes exquisite near a critical point, 
and thermodynamic variables contain a singular depen- 
dence on H and t. The fact that the boundary roughness 
plays the role of external magnetic field is a reflection of 
the fact that small perturbations couple to the flow and 
drive the turbulent state. 

The exponents in the scaling theory given here are 
taken from experiment; however, Gioia and Chakraborty 
have recently pointed out[j^ that the friction factor / 
can be related to the local structure of turbulence, by 
considering the momentum flux at the pipe boundary. 
The Nikuradse data show four features: a hump, the 
Blasius regime, a shallow minimum, and the Strickler 
regime. The scaling argument presented here implies 
that the Blasius and Strickler regimes are both manifes- 
tations of inertial range scaling coupled with wall friction, 
and indeed, Gioia and Chakraborty find from momentum 
flux considerations that this is sufficient to reproduce the 
Blasius and Strickler regimes. If they then include the 
dissipation range in their formula for the friction factor, 
they find that this reproduces the shallow minimum be- 
tween the Blasius and Strickler regime. Interesting their 
formulae for the pure inertial and the inertial + dissi- 
pation ranges both satisfy Eq. The hump in the 
friction factor arises from the energy-containing range, 
lies outside the range bracketed by Blasius and Strickler 
and is absent from their predictions if this spectral fea- 
ture is not included in their formula. Thus, in summary, 
the features which are described by Eq. are indeed 
related to the spectral features of the inertial range (and 
the dissipation range). Thus, just as in critical phenom- 
ena, the large-scale phenomenology of turbulence can be 
related to the power law fluctuations. Presumably, the 
scaling result Eq. Q can be derived from a renormaliza- 
tion group argument 9]. 

Eq. Q describes the nonequilibrium driven steady 
state of fully-developed turbulence, which is to be con- 
trasted with the phenomenon of the transition from 
laminar to turbulent flow. In the case of turbulent 
pipe flow, it is generally accepted that laminar Hagen- 
Poiseuille flow is linearly stable at all Reynolds num- 
bers, and that there is a finite amplitude instability to 
turbulence [2^. liflj. whose amplitude varies pfflj] as Re^ 1 . 
Thus, this transition has the character of a first order 
transition in thermodynamics, and is not related to the 
theory given here. 

Our results establish boundary roughness as a key el- 



ement in a proper theoretical description of turbulence, 
in the same way that a proper understanding of the fer- 
romagnetic critical point at zero field would not be pos- 
sible without taking into account the behaviour for non- 
zero values of H . Our analysis highlights the need for 
a definitive set of experiments to replicate Nikuradse's 
data set, with a view to g reater precision and removing 
sources of uncertainty^ in both the data set and in 
the widely-used Colebrookjj^ semi-empirical fit[13|- I n 
ongoing work we are exploring the effects of roughness 
in two-dimensional turbulence |32(, where predictions for 
the analogues of the Blasius and Strickler laws reveal in- 
teresting differences from the three dimensional results. 

I acknowledge valuable discussions with G.I. Baren- 
blatt, Gustavo Gioia and Pinaki Chakraborty, whom I 
also thank for assistance with the figures. This work was 
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through grant number NSF-EAR-02-2I743. 
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